Measurement of the ^^Ne(a,po)^^Na reaction cross section in the burning energy 

region for X-ray bursts 
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The ^*Ne(Q,p)^^Na reaction provides one of the main HCNO-breakout routes into the rp-process 
in X-ray bursts. The ^®Ne(Q,po)'^^Na reaction cross section has been measured for the first time 
in the Gamow energy region for peak temperatures T~2GK using the time-reversal approach in 
inverse kinematics. The astrophysical rate for ground-state to ground-state transitions was found 
to be a factor of 2 lower than Hauser-Feshbach theoretical predictions. Our reduced rate will affect 
the physical conditions under which breakout from the HCNO cycles occurs via the i**Ne(a,p)2iNa 
reaction. 



PACS numbers: 25.60.-t, 26.30.Ca, 25.45.Hi 



Type I X-ray bursters (XRBs) exhibit brief recurrent 
bursts of intense X-ray emission and represent a frequent 
phenomenon in our Galaxy. Recent observations from 
space-borne X-ray satellites (BeppoSAX, RXTE, Chan- 
dra and XMM Newton) have provided a great wealth of 
data and have marked a new era in X-ray astronomy. 
Yet, to fully exploit these observations, similar advances 
in our understanding of the nuclear reactions responsible 
for the bursts are required. XRBs are driven by a ther- 
monuclear runaway on the surface of a neutron star that 
accretes H- and He-rich material from a less evolved com- 
panion in a semi-detached binary system [l[ . Depending 
on the mass accretion rate, high enough temperatures 
and densities can be achieved that trigger hydrogen burn- 
ing through the hot, /3-limited CNO cycles (HCNO) and 
the subsequent ignition of the triple-a process. However, 
the thermonuclear runaway requires a breakout from the 
HCNO cycle and the ignition of the rapid-proton capture 
process (rp-process) at peak temperatures T~ 1-2 GK. 
The ^*Ne(a,p)^^Na reaction is believed to provide the 
mam breakout route at T>0.8GK and p >IQ^ g/cm^ [2|, 
but the actual physical conditions at which the break- 
out occurs depend critically on the accurate knowledge 
of the "'^*Ne(a,p)^"'^Na reaction rate. A direct investiga- 
tion of this important reaction is severely hampered by 
the low intensity (< lO^pps) of radioactive -'^^Ne beams 
presently available and by the further complications as- 
sociated with the use of a ''He gas target. Thus, the 
only two direct measurements available to date extend 
to minimum energies of Ecni= 2.0 MeV Q and Ecin= 1-7 
MeV Q ■ These are still too high compared to the energy 
region Ecm <1.5 MeV of interest for HCNO breakout in 
X-ray bursts. 

The first theoretical estimates of the ^*Ne(Q;,p)^^Na re- 
action rate Q were based on sparse experimental infor- 
mation on the level structure of the compound nucleus 



^^Mg above the a-particle threshold at 8.14 MeV. Inspec- 
tion of the structure of the mirror nucleus ^^Ne reveals 
a high level density at comparable excitation energies 
and suggests that a statistical approach might provide 
a reliable estimate of the reaction rate. However, only 
natural-parity states in ^^Mg can be populated by the 
^^Ne-|-a channel and the resulting level density may be 
significantly smaller than required for a statistical ap- 
proach. Thus, a range of experimental studies have fo- 
cused on the investigation of a-unbound natural-parity 
states in ^^Mg using reactions such as ^^C(^®0,^He)^^Mg, 
25Mg(3He,6He)22Mg, ^'^Mg[^Yle,^}ief^Mg More 
recently, several a-unbound states were identified using 
the ^^Mg(p,t)^^Mg reaction and the improved precision 
on the measured excitation energies resulted in smaller 
uncertainties on the "'*Ne(Q;,p)^"'Na reaction rate 0. Yet, 
a comparison of all the reaction rates currently available 
shows discrepancies of up to several orders of magni- 
tudes both below and above T '^IGK (see Q and ref- 
erences therein). In addition, it remains unclear whether 
the Hauser-Feshbach statistical calculations provide a re- 
liable estimate of the "'^Ne(a,p)^^Na rate in the whole 
temperature region relevant to nucleosynthesis in X-ray 
bursts. 

In this Letter, we report on the results of a time- 
reversal investigation to determine the ^*Ne(Q;,po)^^Na 
reaction cross section at Ecm ~1.2-2.6 MeV. This is the 
first measurement of this critical reaction in the Gamow 
energy region for a maximum peak temperature T~ 
2GK. As only ground-state to ground-state transitions 
can be accessed from the time-reversal reaction, the in- 
ferred cross section for the direct ^®Ne(Q;,p)^"'Na reaction 
represents a lower limit and further investigations will be 
needed to assess any additional inelastic contributions 
(see ^ for such a study in the case of -'*0(q;,p)^''F). 
It should be noted, furthermore, that given the non-zero 
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spins of the interacting nuclei in the entrance channel 
(J'' = 3/2+ and J'' = 1/2+ for ^^Na and proton, respec- 
tively) both natural- and non-natural-parity states can 
be populated in ^^Mg. However, detection of ^^Ne nu- 
clei in their ground state ensures that only natural-parity 
states in ^^Mg, i.e. those of astrophysical interest, have 
been populated by the inverse reaction. 

The time-reversal measurement was carried out at the 
ISAC II facility (TRIUMF) in inverse kinematics us- 
ing a radioactive ^^Na beam incident on (CH2)„ targets 
(thicknesses 311 and 550 /zg/cm^). Six beam energies 
were explored in the region 4.12 MeV/A < Ebeam <5.48 
MeV/A. Typical beam intensities were limited to a max- 
imum of 10^ pps in order to keep the Rutherford scat- 
tering yield in the most forward detectors (see below) 
at an acceptable rate. Heavy ions and a particles were 
detected in coincidence using an array of Double Sided 
Sihcon Strip Detectors (DSSSDs). Detection of '^He ions 
was primarily achieved with a AE-E telescope consisting 
of two MSL |ll| type S2 detectors (65 /xm and 500 ^m 
thick, respectively, each segmented into 48 annular front 
strips and 16 rear sectors) covering a laboratory angu- 
lar range of 9a = 7.0° — 20.2° (the maximum emission 
angle for a particles was less than 20° at all energies in- 
vestigated). Detection of ^^Ne ions was achieved with 
a AE-E telescope using a "CD" detector consisting of 4 
quadrants (each segmented into 16 front strips and 24 
rear sectors) and a "PAD" detector, consisting of 4 un- 
segmented quadrants [llj. The CD- PAD (35 fiia and 
1500 /im thick, respectively) covered a laboratory angu- 
lar range 9m = 1.6° — 6.6° (the maximum emission angle 
for ^^Ne ions was less than ~ 5° at all energies investi- 
gated) . The coincidence detection efficiency was between 
13% and 27%, depending on beam energy, as determined 
by a Monte Carlo simulation assuming an isotropic dis- 
tribution in the center of mass (l^ . 

The energy calibration of the detectors was performed 
independently for each strip using a mixed three-peak 
a source. In addition, the Rutherford elastic scatter- 
ing of 2iNe off at Eiab=5.357 MeV/A was used 
as a fourth calibration point for the CD quadrants 
only. Good events, corresponding to ''He and -"^^Ne ions 
detected in coincidence, were extracted from the raw 
data by imposing appropriate conditions of two-body 
co-planarity, identification of ''He and '^Ne kinematics 
loci, and total energy reconstruction. An example of 
the measured and simulated a-particle kinematics loci 
for the 2iNa(p,a)^^Ne reaction at Eboam= 5.476 MeV/A 
is shown in Figure [TJ The upper and lower loci corre- 
spond to '^Ne nuclei being left in their ground state or 
first excite state (Ea;=1.89 MeV), respectively. This lat- 
ter channel is open at beam energies 5.476 and 4.910 
MeV/ A (and closed at the other energies investigated). 
No events were observed in the lower locus at either beam 
energy, thus indicating that the inelastic channel is not 
strongly populated in the reaction. Experimental data 
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FIG. 1. Simulated alpha-particle kinematics curves {9 vs. 
E) superimposed on experimental data (black squares) from 
the ^^Na(p,a)^*Ne reaction at Ebeam= 5.476 MeV/A. The 
upper (lower) locus corresponds to reactions forming ^^Ne in 
its ground- (first-excited) state. 



were found to be in excellent agreement with simulations 
at all beam energies. 

The ^'Na(p,a)'®Ne cross section was determined from 
the measured coincident yields under the assumption of 
a thin-target condition at each investigated energy, using 
the following equation: 



CT(^cff) 



Y 



where E^g is the effective beam energy at mid-target, Y 
is the measured ^'Na(p,a)'*Ne coincident yield, Nf, the 
total number of incident particles, Nt the number of 'H 
nuclei per unit area in the (CH2)„ target, ^ the coincident 
detection efficiency, and r is the data acquisition (DAQ) 
live time (typically 65%). The number of incident par- 
ticles was calculated from the yield of Rutherford scat- 
tered ^'Na ions off ^^C in the target as detected in the 
innermost strips of the CD detectors, i.e well within the 
grazing angle at each beam energy. 

A summary of experimental yields and cross sections 
for each of the six beam energies investigated is given 
in Table ID Note that the errors in the effective inter- 
action energy represent half the target thickness (in the 
center-of-mass system) at each beam energy. The er- 
rors in the yields and cross sections are statistical only 
and were calculated using Poisson statistics for the two 
highest yields and using the Feldman-Cousins method 
|13| for low statistics for all other yields (68% confidence 
level and zero background assumption). The estimated 
systematic uncertainty of the cross section is ~16% and 
is dominated by uncertainties in the number of projec- 
tile nuclei (8% at all beam energies) and in the number of 
target nuclei (about 8%, mostly due to the uncertainty in 
the energy loss calculations by SRIM2008 [lj|). Changes 
to the detection efficiency due to non-isotropic distribu- 
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TABLE I. Summary of experimental results for each of the six beam energies investigated. 



Ebcam 

(MeV/ A) 


Ecm (P,a) 

(MeV) 


Ecm(a,P) 

(MeV) 


He+^Ne 
Yield 


C 
(%) 


Nb'' 
particles 


(T(p,a) 
(mb) 




a(a,p) 
(mb) 


5.476'' 


5.21±0.06 


2.57±0.06 


33±6 


19 


(2.9 ±0.3) X 10"' 


0.35±0.06 




1.7±0.3 


4.910'' 


4.61±0.12 


1.97±0.12 


0+3.3 
°-2.7 


27 


(3.2 ±0.3) X 10^" 


(3.0±1;?,) X 10- 


2 


17+007 


4.642^ 


4.40±0.07 


1.76±0.07 


23±5 


27 


(7.1 ±0.7) X 10" 


(5.3 ± 1.1) X 10 


-3 


(3.1 ±0.6) X 10-2 


4.619*^ 


4.32±0.12 


1.68±0.12 




25 


(5.0 ±0.5) X 10" 


(3.8ljj) X 10- 


3 


(2.3«:I) X 10-2 


4.310'' 


4.02±0.13 


1.38±0.13 


4+2.8 
^-1.7 


16 


(1.47 ±0.15) X 10^2 


(5.6t^-^) X 10- 


4 


(3.8+2-^) X 10-3 


4.120'' 


3.83±0.13 


1.19±0.13 


n+2.3 
^-1.3 


13 


(6.9 ±0.7) X 10^2 


{7A+_il) X 10- 


5 


(5.5+«;2) X 10-^ 



^ 311^ig/cm2 (CH2)n target 
'' 550^ig/cm2 (CH2)„ target 

as determined by Monte Carlo simulation 
'' as determined by Rutherford elastic scattering of 2^Na beam off ^^C nuclei in the target (see text). 



tions (up to ^=3) were also explored by Monte Carlo sim- 
ulation. Deviations from the isotropic case amount to at 
most 20% at the highest beam energies and to at most 
50-60% at the two lowest beam energies (depending on 
the actual I value) and are therefore comparable to, or 
smaller than, the quoted statistical uncertainties. 

The ^*Ne(Q;,po)^^Na reaction cross section was inferred 
by using the principle of detailed balance [l5| and is 
shown in Figure [2] as a function of (q;,p) center-of-mass 
energy. The figure also shows the theoretical predictions 
based on Hauser-Feshbach calculations [l^ for ground- 
state to ground-state transitions (hereafter HFgs) and 
ground-state-to-all-states transitions (HFaii). Since our 
data only provide a lower limit to the cross section, 
the comparison with theoretical predictions is made only 
with HFgs. Surprisingly, a good agreement is found for 
the two lowest energy points, while a discrepancy of up to 
a factor of 2 is observed at the highest measured energy. 
This is contrary to expectations, as lower energies cor- 
respond to lower excitation energies, and therefore lower 
level densities, in the compound nucleus. The reason for 
such a trend is at present not understood. According to 
HF predictions, contributions to the total cross section 
from the inelastic channels in the energy region investi- 
gated amount to at most '^3x those for the ground-state 
to ground-state transition. 

The astrophysical ^^Ne(Q;,po)^^Na reaction rate as a 
function of temperature was calculated by numerical in- 
tegration of our experimental cross section data and is 
shown in Figure [3l A comparison with the rate of the 
direct measurement of Groombridge et al. reveals 
significant discrepancies both in energy dependence and 
magnitude, with a disagreement of up to a factor of ~ 25 
at T=2.4GK. It should be pointed out that such discrep- 
ancies would persist (albeit to a lesser extent) even by 
allowing for inelastic contributions to our rate according 
to HF predictions. The origin of this large discrepancy is 
not clear but we note that the resonance widths reported 
in 0, 17 1 imply unphysically large alpha-decay partial 
widths. Our rate is also between a factor of about 10-40 
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FIG. 2. Experimental ^*Ne(a,po)2^Na reaction cross section 
(black dots) as a function of Ecm''' . Predictions based on the 
Hauser-Feshbach calculations [3 for ground-state to ground- 
state transitions (full line) and ground-state to all-states tran- 
sitions (dashed line) are also shown for comparison. 



lower than the rate deduced using an indirect approach 
by Matic et al. 9] for 0.9 GK<T<2.4 GK. Such indi- 
rect methods, while a valuable guide, are subject to large 
uncertainties in the choice of resonance strength values, 
which are not determined experimentally. Our reaction 
rate agrees well with the overall energy dependence of 
the HFgs rate but is typically a factor of '^2-3 lower in 
the whole temperature region (T=1.0-2.4 GK) that cor- 
responds to our measured energy range. Outside of this 
temperature range our rate can be extrapolated but its 
reliability decreases drastically, especially for T < 0.95 
GK and thus a direct comparison with the rate of 6] is 
not possible. 

Based on our results, we can infer that the breakout 
from the HCNO cycle via the ^^Ne(a,p)^^Na reaction is 
delayed and occurs at higher temperatures than previ- 
ously predicted. However, more detailed conclusions re- 
quire full hydro-dynamical model calculations for X-ray 
bursts that are beyond the scope of this Letter. 
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date and, for the first time, within the Gamow energy re- 
gion of this reaction in X-ray bursts. Our results indicate 
that a breakout from the HCNO via the ^8Ne(a,p)2iNa 
reaction should occur at higher temperatures than pre- 
viously assumed. 
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FIG. 3. (Colour online) Astrophysical ^*Ne(a,p)^^Na reac- 
tion rate from the present work (solid black line) calculated 
by numerical interpolation of our cross-section data. Rates 
from (black dotted line) and ^ (blue dash-dotted line) 
are larger than the present rate by up to a factor of 25 and 
40, respectively at T=2.4GK. Theoretical rates based on HP 
formalism [l^ (red dashed line: ground-state to ground-state 
transitions; dotted red line: ground-state to all-states transi- 
tions) are also shown for comparison. 



TABLE II. Calculated rates of the ^*Ne(a,p)^^Na reaction as 
a function of temperature. 



Temperature Na 
Tg (K) limit low 


< av > (cm^ 
limit high 


mor^ s-1) 
arithmetic mean 


0.95 


8.4x10 


-4 


3.1x10" 


-3 


(2.0 


± l.l)xlO" 


-3 


1.05 


4.3x10 


-3 


1.5x10" 


-2 


(9.8 


± 5.5)xl0" 


-3 


1.15 


1.7x10 


-2 


5.8x10" 


-2 


(3.8 


± 2.0) X 10" 


-2 


1.25 


5.7x10 


-2 


1.8x10" 


-1 


(1.2 


± 0.6) X 10" 


-1 


1.35 


1.6x10 


-1 


4.9x10" 


-1 


(3.3 


± 1.7)xl0" 


-1 


1.45 


4.1x10 


-1 


1.2 




(7.9 


± 3.9) X 10" 


-1 


1.55 


9.2x10 


-1 


2.6 






1.8 ± 0.8 




1.65 


1.9 




5.2 






3.6 ± 1.6 




1.75 


3.8 




9.8 






6.8 ± 3.0 




1.85 


6.9 




1.7x10^ 


-1 


(1.2 


± 0.5) xlO^ 


-1 


1.95 


1.2x10 


f 1 


3.0x10"' 


-1 


(2.1 


± 0.9) xlO^ 


-1 


2.05 


2.0x10 


f 1 


4.8x10^ 


-1 


(3.4 


± 1.4) xlO^ 


~1 


2.15 


3.1x10 


f 1 


7.4x10^ 


-1 


(5.3 


± 2.2) xlO^ 


-1 


2.25 


4.7x10 


f 1 


1.1x10"^ 


'2 


(7.9 


± 3.2) xlO"' 


-1 


2.35 


6.9x10 


f 1 


l.exlO"' 


-2 


(1.2 


± 0.5) X 10+^ 


2.45 


9.8x10 


f 1 


2.3x10^ 


-2 


(1.6 


± 0.7) xlO^ 


-2 
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